• GRUMMAN AEROSPACE CORPORATION 
LM MEMORANDUM 


From: R. Schindwolf [i. 

To: Distribution 

Subject: Rate of Descent (ROD) Mode Stability and Performance Analysis. 
SUMMARY': 

The results of this analysis indicates that the stability of the ROD 
control system can be optimized by properly matching the LGC programmed 
gain constants to the system time delays. For Apollo 12 the engine response 
time assumed in choosing the gain.constants was longer than the actual 
engine response time. This resulted in a lightly damped system and in con¬ 
junction with the extraneous accelerations introduced by the IMU mounted 
off the C.G. resulted in the large throttle oscillations observed during 
the Apollo 12 landing. The throttle oscillations observed during FCI lab. 
tests were due to accelerometer noise introduced by quantization in the 
PIPA math model. To reduce the sensitivity of the loop to extraneous accel¬ 
erations it is recommended that the ICC gain constants and \be changed 
to 0.2 and .1 seconds respectively. These gains will result in a proper match 
to the best estimate of the system time delays. In addition it is recommended 
that the program change suggested by MIT.which approximately compensates for 
the extraneous accelerations introduced by the off-C.G. IMU be incorporated. 
These changes will reduce the peak to peak throttle variations from 2500 pounds 
as observed on Apollo 12 to approximately 300 - 400 pounds. 

INTRODUCTION : 

. FMES/FCI powered descent tests indicated that throttle oscillations in the 
order of 800 pounds peak to peak occurred during the ROD mode. A review of the 
Apollo 12 flight data indicated that throttle oscillations as large as 2500 
pounds peak to peak occurred during ROD operation. The purpose of this analysis 
was to determine if these oscillations were due to an inherent instability in 
the ROD control loop or to extraneous forcing functions. 

DISCUSSION : 

A loop diagram of the ROD throttle loop is shown in Figure 1. This loop 
was derived from the Guidance Equation Section of the GSOP. The e” ST term 
used in the transfer functions represents a time delay of one computation 
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period (l second). Due to the sample data nature of the system it is more 
convenient to use Z transforms instead of Laplace transforms (Z - e ST ). The 
loop diagram using Z transforms is shown in Figure 2. 

The Z transform from A A* to H is not readily apparent but can be deduced 
by physical reasoning. The Z transform of a linear system is simply the 
response of the system at the sampling times to an input impulse. A unit 
impulse in A A* will result in an H ramp of unit slope and delayed by T^ 



The Z transform is therefore 
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System stability is determined by the location of the roots of 1 + G (z) 
where G (Z) is the open loop gain function . For the system to be 

stable the roots of 1 + G (z) must lie inside the unit circle in the Z plane. 
G (Z) can be obtained from the loop diagram and is as follows: 
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where T. 


D 


total system time delay including computational 
time and engine delays. 


T 1 
1 D 


Computer estimate of total time delay 
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Therefore 

I f [ 7 - + Qi- T o'] + (l + Hr)^ ') i 1 ' 7 *) 2 + 
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and system stability is determined by the zeros of the polynomial 


(z-lf [i t(l* ij-K] +('/(- J0(i-l)§-%)^ f r,] t Jr z[(i-T,)2*T,] 

Rearranging terms this polynomial can be written as follows: 

2 a :-fz-^ L )+ J^z2°.(z-i)2 +(To-7 D Yi + ^)(2-\f 

J ^ 

If LAG and Tj* are properly matched to the system delays (LAG = = T ) the 

last two terms drop out and the characteristic equation is 

2 Vz-^) 

With CT =1.5 this becomes 

2 f^a-. 22 ) 

Obviously this system is extremely stable with a double root at the origin 
and a root at + . 33 . 


In order to determine the effects of mismatches the root loci of Figures 
3, b and 5 were drawn. Figure 3 gives the root locus as LAG is varied from 
+1.2, In this diagram it is assumed that the system time delay y T , 
is .3 seconds and that is properly matched. From this locus it can be 
seen that any value of LAG other than .3 (the matched value) will result in 
the root moving away from the origin toward the'unit circle and will be less 
stable. ^ Figures b and 5 give the same loci except with mismatched by + .1 
seconds (over compensation) and -.1 seconds (under compensation) respectively. 
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It can be seen that the T^ mismatch (both positive and negative) also 

causes a reduction in stability since the roots are further from the 

origin than in Figure 3» It is therefore concluded that by choosing 

LAG and T-p to be equal to the best estimate of the system time delay, 

optimum stability will be achieved. Since the best estimate of total 

system time delay is approximately .3 seconds (.2 seconds computation 

tine and- .1 seconds throttle lag) the — te - rm in the pgc p rogr am 

should be 0.2. The LGC determines T^ by measuring the computational 

delay and adding a prestored value of throttle lag (O',,). For optimum 

tn 

stability this stored value should equal the best estimate of throttle 
lag (.1 seconds). 

LAG 

The stored values of ~ and O' th used on Apollo 12 were 0.4 and 
0.2 seconds respectively. These values result in the root locations 
circled in Figure 4. It can be seen that one of the roots is close to 

It is therefore 

be changed 


'V 
v th 


the unit circle resulting in a lightly damped system, 
recommended that for future flights the value of and 

to 0.2 and 0.1 seconds respectively. 

In order to explain the large throttle variations the response of the 
system to accelerometer noise was computed. The closed loop transfer 
function relating acceleration command to accelerometer AV errors is as 

f01 T2 fr-.) $+ & **) * - (1+ i^)l 

The transfer functions using the properly matched gains and the Apollo 12 
gains are as follows: 
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The acceleration command resulting from a one= : foot per second velocity 
error impulse is obtained by expanding these transfer functions as power 
series ing- . The two responses are plotted in Figure 6. 

The Apollo 12 response is particularly significant since the magnitude: 
of the velocity errors introduced by the.off-CG IMU was in the order of 
1 ft/sec. From the plot it can be seen that the peak to peak acceleration 
command variation is as large as 5*9 ft/sec 2 -which is equivalent to 3000 lbs 
of thrust variation. This agrees quite well with the thrust variations 
observed on Apollo 12. In addition, the respoilse indicates poor damping as 
expected from the stability analysis. The response for the properly matched 
system shows similar peak amplitudes but decays rapidly after the first two 
sample periods. 

These response characteristics can also be used to explain the throttle 
variations observed in FMES/FCI lab tests. In the FMES there are no velocitv 
errors due to an off-C.G. IMU but errors are introduced by 50 millisecond time 
quantization of the simulated accelerometer outputs. This can result in 
instantaneous velocity errors as large as .25 ft/sec and peak to peak throttle 
variations in the order of 700 v pounds for the Apollo 12 gains. This agrees 
quite well with the results observed in tests* 

The throttle variations in the ROD mode can be reduced by either com¬ 
pensating for the erroneous accelerometer data or by changing the control 
loop to make it less, sensitive to short term acceleration errors. A program 
modification that compensates for the errors due to the off-C.G. IMU has been 
developed by MIT and tested in the FCI lab. The test results indicated throttle 
variations of kOO lbs under conditions similar to those experienced in Apollo 12. 
In order to reduce the system sensitivity it is desirable to make the gain changes 
mentioned previously to match the system time delays. This does not change the 
initial amplitude appreciably but it does reduce the settling time. Other 
values of LAG and throttle lag time may result in somewhat smaller initial 
response but will result in reduced damping and are therefore not desirable. 

One method of reducing the initial amplitude and still maintaining good damping 
is to increase CT . The response to a unit velocity impulse for 0" = 3.0 is 
shown in Figure 7« It can be seen that the peak to peak response is reduced 
by about 30$. The disadvantage of this gain change, however, is that it 
results in increased response time to ROD inputs. 
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In order to obtain a significant reduction in accelerometer noise 

sensitivity without increasing ROD response time it is necessary to 

modify the basic loop configuration. The most significant feedback 

paths for erroneous acceleration data is through the term (lAFj + j£l&) 

' 1 / 

vhich is used in the P66 guidance routine and in the throttle command 
routine.. This term is essentially the present value of vehicle thrust 
acceleration. Rather than using accelerometer derived data this term 
can be replaced by the commanded acceleration from the previous cycle. 

. All other gains and parameters would be unchanged. This will eliminate 
two accelerometer feedback paths ano therefore reduce considerably the 
response to accelerometer noise. The response to ROD inputs and system 
stability will not be affected by this change. The response of this 
configuration to a one ft/sec accelerometer error is shown in Figure 8 
and it can be seen that the peak to peak response has been reduced by a 
factor of about 4. This is a significant reduction and it is recommended 
that the program changes necessary to implement this technique be evaluated 

C * . 

CONCLUSIONS : 

The throttle oscillations observed on Apollo 12 were due to the control 

system beifrg lightly damped and being forced by exteraneous accelerations 

due to the off~C.G. IMU. To improve the. damping and stability of the system 

it is recommended that the LGC program constants and O' , be changed 

th 

to .2 and .1 respectively. In addition, the MIT program change which auprox- 
imately compensates for the offset IMU accelerations should be incorporated. 
The alternate loop configuration suggested should be studied further since it 
the loop response to any other extraneous acceleration data such 
as PIPA quantization, time quantization, vibration, etc. 

RS/sjd 

cc: L. Tucker 
‘ C. Tillman 
J. Marino Sy\ 

M. Rimer « 

P. Hoffman 
V/. Holmer 
\ S. Greene 

P. Kelly 

C.~Hackler MSC/GCD 
R. Steele" — 

W. Nufer 
. G. Sullivan 
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FIGURE 4 - Z PLANE HOOT LOCUS (T. 
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FIGURE 5 - Z PLANE ROOT LOCUS (T* - T = -.1) 
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FIGURE 7 .. RESPONSE TO 1 FT/SEC ACCELEROMETER IMPULSE 
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FIGURE 8 - RESPONSE TO 1 FT/SEC ACCEIEROMETER IMPULSE 
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